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Measurement of the Branching Fraction of Meson Decay to 

a^(1260) 7f- 



Abstract 

We present a preliminary measurement of the branching fraction of the B meson decay 

a^"(1260)7r- with 0^1^(1260) vr+vr+vr-. The data sample corresponds to 218 x 10^ BB 
pairs produced in e~^e~ annihilation through the T{4S) resonance. We find the branching fraction 
(40.2 ±3.9 ±3.9) X 10"^, where the first error quoted is statistical and the second is systematic. The 
fitted values of the ai(1260) parameters are ruai = 1-22 ± 0.02 GeV/c^ and Ta^ = 0.423 ± 0.050 
GeV/c2. 
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1 INTRODUCTION 



We report on the preliminary measurement of the branching fraction a^(1260)7r~ with 

a^(1260) 7r^7r"^-7r~ jJI • The ai(1260) — > Svr decay proceeds mainly through the intermediate 
states (7r7r)p7r and (vrvrjo-TT [S] • 

The study of this decay mode is complicated by open questions on the parameters of the 
ai(1260) meson. There are large discrepancies between these parameters when comparing results 
from analyses involving hadronic interactions jSI and r decays jl]. Therefore, it is important 
to verify the theoretical prediction of the branching fraction for this decay mode and have new 
measurements of the ai(1260) parameters. A theoretical calculation of the branching fraction of 
this decay mode has been made by Bauer, Stech and Wirbel (BSW) ^ within the framework of the 
factorisation model. They predict a value of 38 x 10^ , assuming -vf^ = 0.08. It is also important 

to note that the a^(1260)7r~ channel can be used to measure the Cabibbo-Kobayashi- 

Maskawa angle a of the Unitarity triangle |^ . We presented a preliminary version of this analysis 
at ICHEP'04 7 , using an integrated luminosity of 112 fb~^ and the measured branching fraction 
was (42.6 lb 4.2 it 4.1)xl0~®. For the branching fraction of B^ a^(1260)7r~ an upper limit of 
49 X 10~^ at the 90% confidence level (C.L.) has been set by CLEO collaboration 8^ while the 
DELPHI collaboration has set the 90% C.L. upper limit of 28 x 10~^ for the branching fraction 
of B^ 47r. 

Below we present the details of the analysis for the measurement of the branching fraction for 
B^ a^(1260)7r~ — > 27r"'"27r~. Presently, we do not distinguish between the final states (7r7r)p7r 
and {Tnr)aTT. Such an analysis would require a study of the angular distributions of the decay 
products. Possible background contributions from B^ decays to a2^(1320)7r^ and 7r^(1300)7r~ are 
studied and taken into account while in the preliminary version presented at ICIIEP'04 they were 
neglected. 



2 THE BABAR DETECTOR AND DATASET 

The results presented in this paper are based on data collected in 1999-2004 with the BABAR 
detector ^JJI at the PEP-II asymmetric e~^e~ collider jT^ located at the Stanford Linear Accelerator 
Center. An integrated luminosity of 198 fb~^, corresponding to 218 million BB pairs, was recorded 
at the T{4:S) resonance ( "on-resonance" , center-of-mass energy ^/s = 10.58 GeV). An additional 
15 fb^^ were taken about 40 MeV below this energy ("off-resonance") for the study of continuum 
background in which a light or charm quark pair is produced instead of an T{4S). 

The asymmetric beam configuration in the laboratory frame provides a boost of /?7 = 0.56 to 
the T(4S). Charged particles are detected and their momenta measured by the combination of 
a silicon vertex tracker, consisting of five layers of double-sided silicon microstrip detectors, and 
a 40-layer central drift chamber, both operating in the 1.5-T magnetic field of a solenoid. The 
tracking system covers 92% of the solid angle in the center-of-mass frame. 

Charged-particle identification is provided by the average energy loss (dE/dx) in the tracking 
devices and by an internally reflecting ring-imaging Cherenkov detector (DIRC) covering the central 
region. A K/ir separation of better than four standard deviations (a) is achieved for momenta below 

3 GeV/c, decreasing to 2.5 a at the highest momenta in the B decay final states. Photons and 
electrons are detected by a CsI(Tl) electromagnetic calorimeter while muons are identified in the 
magnetic flux return system. 
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3 ANALYSIS METHOD 



Monte Carlo (MC) simulations jT2| of the signal decay mode, of continuum and BB backgrounds 
are used to establish the event selection criteria. We make several particle identification re- 
quirements to ensure the identity of all signal pions. For the bachelor charged track we re- 
quire an associated DIRC Cherenkov angle between —2a and +5cr from the expected value 
for a pion. A B meson candidate is characterized kinematically by the energy-substituted mass 
TTiES = \J s + Po ■ VbY /Eq — and energy difference A£' = E*^ — \\/s, where the subscripts 
and B refer to the initial T{AS) and to the B candidate in the lab-frame, respectively, and the 
asterisk denotes the T{AS) frame. We require \/\E\ < 0.2 GeV and 5.25 < tues < 5.29 GeV/c^. We 
select a||^(1260) candidates with the following requirement on the invariant mass: 0.8 < < 1.8 
GeV/c^. The intermediate dipion state is required to have an invariant mass between 0.51 and 1.1 
GeV/c^. The momentum of a|,^(1260) in the center-of-mass frame is required to be between 2.3 
and 2.7 GeV/c. To reduce fake B meson candidates we require p(x^) > 0.01 for the B vertex fit. 
The angular variable Tia^ (cosine of the angle between the direction of the bachelor vr and the flight 
direction of the B in the ai(1260) meson rest frame) is required to be between —0.85 and 0.85 to 
suppress combinatorics. 

To reject continuum background, we make use of the angle 6t between the thrust axis of the 
B candidate and that of the rest of the tracks and neutral clusters in the event, calculated in 
the center-of-mass frame. The distribution of cosOt is sharply peaked near ±1 for combinations 
drawn from jet-like qq pairs and is nearly uniform for the isotropic B meson decays; we require 
|cos0t| < 0.65. The remaining continuum background is modelled from "off-resonance" data. 
We use Monte Carlo simulations of B^B^ and B^B~ decays to look for BB backgrounds, which 
can come from both charmless and charm decays. We find that the decay mode B^ — > D~'ir~^, 
with D~ —>■ K^'K~'K~ and D~ — > K^'k~ , is the only significant background. It is included in the 
maximum likelihood fit. Final results have been corrected for a small background contribution due 
to charmless decays. 

We use an unbinned multivariate maximum-likelihood fit to extract the signal yields for 
S° at{l2m)TT-. The likelihood function incorporates five variables. We describe the B decay 
kinematics using: Ai?, msS) ^a^^ a Fisher discriminant J- , and an angular variable A. The Fisher 
discriminant combines four variables: the angles in the T{AS) frame of the B momentum and B 
thrust axis with respect to the beam axis, and the zeroth and second angular moments -Lo,2 of the 
energy flow around the B thrust axis. The moments are deflned by 

Lj = ^Pi\cosei\^ , (1) 

i 

where pi is the momentum of particle i, 9i is the angle between the direction of particle i and the 
trust axis of the B candidate and the sum excludes tracks and clusters used to build the B candidate. 
We have used an angular variable A in order to distinguish o^(1260)7r~ from a2^(1320)7r~ and 
7r+(1300)7r-. If X is our resonance ai(J^ = 1+), a2(J^ = 2+) or 7r(1300)(J^ = 0") that decays 
into three pions, we evaluate in the X meson rest frame the cosine of the angle between the normal 
to the plane of the three pions and the flight direction of the bachelor pion. Since we have on 
average 1.5 B candidates per event, we choose the best one using a quantity computed with 
the p mass. Since the maximum correlation between the observables in the selected data is 4%, 
we take the probability density function (PDF) for each event to be a product of the PDFs for 
the separate observables. The product PDF for event i and hypothesis j, where j can be signal (3 
types), continuum background or BB background, is given by 
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V] = VjimEs) ■ r,{AE) ■ Vj{T) ■ Vj{ma,) ■ Vj{A). (2) 

There is the possibility that a track from a signal candidate is exchanged with a track from 
the rest of the event. We call these events "self-cross-feed" (SCF) events. The fraction of SCF 
events with respect to the total number of signal events for each type k of signal, fscF^j is fixed 
to the value found with Monte Carlo signal events (26%). The likelihood function for the event i 
is defined as : 

3 

k=l 

where nk{k = 1,3) is the yield for a||''(1260) vr^, a^(1320)7r^, and 7r+(1300) 7r^ respectively, Uqq 
the number of continuum background events, n^^^ the number of BB background events D~7r^ 
with D~ — > K~^7:~7r~ and ?t-b^2 number of BB background events D^tt^ with D~ Kg7r~. 
The extended likelihood function for all events is : 

« i 

where rij is the yield of events of hypothesis j found by the fitter, and N is the number of events 
in the sample. The first factor takes into account the Poisson fluctuations in the total number of 
events. 

We determine the PDFs for signal and BB backgrounds from MC distributions in each observ- 
able. For the continuum background we establish the functional forms and initial parameter values 
of the PDFs with off-resonance data. We allow the signal oi(1260) PDF parameters and the most 
important qq background PDF parameters to float in the final fit. The distributions of invariant 
mass of ai(1260), 02 (1320) and 7r(1300) in signal events are parameterized as relativistic Breit- 
Wigner line-shapes with a mass dependent width which takes into account the effect of the mass 
dependent p width. The rriES ^-nd AE distributions for signal are parameterized as double gaussian 
functions. Slowly varying distributions are parameterized by linear functions. The combinatoric 
background in mss is described by a phase-space- motivated empirical function J^]- We model the 

distribution using a Gaussian function with different widths above and below the mean. The 
A distributions are modelled using Gaussians in (1260) vr" and polynomials in 02^(1320) vr" and 
vr+(1300)vr-. 

4 RESULTS 

We present the measurement of the branching fraction of the B decay to a^(1260) vr~, considering 
0^(1320) vr~ and vr"'"(1300) vr~ as sources of background. By generating and fitting simulated sam- 
ples of signal and background events, we verify that our fitting procedure is working properly. We 
find that the minimum ln£ value for the on-resonance data lies well within the ln£ distribution 
from these simulated samples. Fits to data show no evidence of vr+(1300) vr~ , since a negative yield 
is obtained for this resonance. For this reason the vr"^(1300) vr^ component has been left out in 
final fits to the yields. 

The reconstruction efficiency is obtained from the fraction of signal MC events passing the selec- 
tion criteria once corrected for a bias detected in the fit yield. This bias (about 6%) is determined 
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Quantity 


a^(1260)7r" 


Signal yield 


867 ± 85 


Reconst. e (%) 


19.8 


US^ (%) 


50 


Stat. sign, (a) 


18.5 


i3(xl0-6) 


40.2 ±3.9 ±3.9 



Table 1: Final fit results in — > 0^^(1260) vr^. Fitted signal yield, the final reconstruction 
efficiency (e), the daughter branching fraction product, the statistical significance, and the central 
value of the branching fraction with statistical and systematic errors. 

from fits to simulated samples, each equal in size to the data and containing a known number of 
signal MC events combined with events generated from the background PDFs. 

The fitted values of the ai(1260) parameters are: = 1.22 ± 0.02 GeV/c^ and Tqi = 

0.423 ± 0.050 GeV/c^ . In Table 1 we show the results of the fits for on-resonance data. The 
statistical error on the number of events is taken to be the change in the central value when the 
quantity — 21n£ changes by one unit. The statistical significance is taken as the square root of 
the difference between the value of — 21n£ for zero signal and the value at its minimum. In Fig. 
^we show the m-ESj Ai?, and A projections made by selecting events with a signal likelihood 
(computed without the variable shown in the figure) exceeding a threshold that optimizes the 
expected sensitivity. 

5 SYSTEMATIC STUDIES 

Most of the systematic errors on the yields that arise from uncertainties in the values of the PDF 
parameters have already been incorporated into the overall statistical error, since they are fioated 
in the fit. We determine the sensitivity to the other parameters of the signal PDF components by 
varying these within their uncertainties. The result is shown in the first row of Table 2. This is 
the only systematic error on the fit yield; the other systematics apply to either the efficiency or the 
number of BB pairs in the data sample. 

The uncertainty in our knowledge of the efficiency is found to be 0.8A'^t%, where Nt is the 
number of signal tracks. We estimate the uncertainty in the number of BB pairs to be 1.1%. 
The fitting algorithm introduces a systematic bias of 2.8%, which was found from fits to simulated 
samples with varying background populations. Published world averages j2j provide the B daughter 
branching fraction uncertainties. The systematic error from ai{1260)K cross-feed background is 
estimated to be 1.4%, while the systematic error due to SCF is found to be 3.5%. We also take 
into account systematic differences between data and MC for the cos^t selection (1.8%) and the 
possibility of interference between the oi and 02 amplitudes (4%). The values for each of these 
contributions are given in Table 2. 
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Figure 1: Projections of mEs(a)) ^E{h), ai mass(c) and A(d) for a^(1260)7r~. Points with errors 
represent data, dotted lines the background from continuum and BB combined, solid curves the 
full fit functions. These plots are made with a cut on the signal likelihood and thus do not show 
all events in the data sample. 
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Quantity 


ajir 


Fit yield 


6.2 


Fit eff/bias 


2.8 


Track multiplicity 


1.0 


Tracking eff 


3.2 


Number BB 


1.1 


SCF 


3.5 


til -A. CIObb-IccQ 


1 A 


MC statistics 


0.6 


COS 6t 


1.8 


ai-a2 Interf. 


4.0 


Total 


9.6 



Table 2: Estimates of the systematic errors (in percent). 



6 SUMMARY 

We have obtained a preliminary measurement of the branching fraction for B^ meson decays to 
a]'"(1260) 7r~ with a5j^(1260) — tt+tt+tt^. The measured branching fraction is: 

B{B^ o|(1260) vr") = (40.2 ± 3.9 ± 3.9) x lO^^ (5) 

The fitted values of the ai (1260) parameters are: nia^ = 1.22±0.02 GeV/c^ and = 0.423±0.050 
GeV/c2. 
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